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Abstract. We establish necessary and sufficient conditions for the existence of 
factorizable steady states of the Generalized Zero Range Process. This process allows 
transitions from a site i to a site i + q involving multiple particles with rates depending 
on the content of the site i, the direction q of movement, and the number of particles 
moving. We also show the sufficiency of a similar condition for the continuous time 
Mass Transport Process, where the mass at each site and the amount transferred in 
each transition are continuous variables; we conjecture that this is also a necessary 
condition. 



PACS numbers: 0250.Ey, 0540.-a 
1. Introduction 

The classical zero range process (ZRP) is a widely studied lattice model with stochastic 
time evolution pQ. To define the process consider a cubic box A C 1< d with periodic 
boundary conditions, i.e. a <i-dimensional torus. At each site i of A there is a random 
integer- valued variable rii G {0, 1, . . .}, representing the number of particles at site i. 
The time evolution is specified by a function a q {rii) giving the rate at which a particle 
from a site i containing n, particles jumps to the site i + q, where q runs over a set 
of neighbors E (the most common choice being E = {±e 1; ±e2, •••})• The name zero 
range indicates the fact that the jump rate from i to % + q depends only on the number 
of particles at i. 

It is easy to see, when the system is finite, that when a q (n) + a_ g (n) > 5 > 
for all n > (we always have a(0) = 0) and for many choices of E (in particular any 
that spans Z d ) then all configurations with a given total particle number N = ^2 ieA U{ 
are mutually accessible, and hence there is a unique stationary measure Pa(jv,N) for 
each N. Normalized superpositions of these measures yield all of the stationary states 
of this system. Conversely, given a measure for which there is a nonzero probability 
of N particles being present in the system one can obtain P\(n; N) by restricting that 
measure to configurations with N particles. 
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The ZRP was first introduced in |2j. It was assumed there and in most subsequent 



works that the rates a q (n) are of the form 



u q {n) = g q a{n) (1) 

with g q independent of n and a(n) independent of q. In this case the system has the 
unique steady state given by E] 




P A (n; AT) = C N 5 Vn 8 - N TTKO ( 2 ) 



where 



p(n) = mSw (3) 



Cn and c are normalization constants given by 

^ A" 

^I1Ll«(*) 



00 \n 



c »=[ e n^) ( 5 ) 

The unique stationary measure Pa{r; N) is thus a restriction to configurations with A" 
particles of the product measure YlieAP( n i) with single-site distribution p(n). 

In the limit A — > Z d with AT/|A| — > p the only stationary extremal measures, i.e. the 
only stationary measures with a decay of correlations, are product measures with p(n) 
given in Equation (j3J) as the distribution of single site ocupation numbers j3 r These 
states are parameterized by A, which plays the role of the fugacity in an equilibrium 
system, with different values of A corresponding to different expected particle densities 
p, where 

n=l n=l - 1 ±K ~ 1 v ' 

Recently there has been a revival of interest in the ZRP. For certain choices of 
a(n), for example when a(n) ~ b/n 2 for large n, the ZRP on Z exhibits a transition 
between a phase where all sites almost certainly contain finite numbers of particles to a 
'condensed' phase where there is a single site containing an infinite number of particles 
[31 This condensation has attracted attention as a representative of an interesting 
class of phase transitions in one dimensional non-equilibrium systems, and has also been 
applied to models of growing networks [HEj- 

Evans, Majumdar and Zia |7j have proposed a generalization of the ZRP, called 
a Mass Transport Model (MTM). They considered a one dimensional lattice on which 
there is a continuous 'mass' rrii > at each site, with a parallel update scheme in which 
at each time step a random mass fa, < \ii < mi moves from each site i e Z to the 
neighboring site i + 1 with a probability density (j){fa mi). This process shares many 
features with the (totally asymmetric) Zero Range Process, in particular the existence 
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of a condensation transition in certain cases [Sj. One very significant difference from 
the ZRP, however, is that the system has a product measure steady state if and only 
if </>(//, m) satisfies a certain condition. This gives rise to a similar condition on the 
sequential update case (which we will call a Mass Transport Process or MTP), where 
the rate of transitions from site i to i+ 1 is given by cn(jU, m). This process has a product 
measure steady state if and only if there exist functions g and p such that 



An interesting question, then, is whether similar criteria for the existence of a 
product measure exist for such processes in higher dimension and with movement in 
both directions allowed. This question is already relevant for the ZRP. A particular 
case in d = 2, studied in [H| HH). has 



i.e. a constant rate (independent of n) per occupied site for moving in the ±x direction 
and a rate proportional to n in the ±y direction. A treatment of this system based 
on fluctuating hydrodynamics and computer simulations (originally conducted on a 
similar but not quite equivalent system, but which we have reproduced on this system) 
suggests that this particular system has correlations between occupation numbers at 
different sites a distance D apart decaying according to a dipole power law D~ 2 . This 
behavior, which is very different from a product measure steady state or its projection 
(J2J), is conjectured to be generic for nonequilibrium stationary states of systems with 
non-equilibrium particle conserving dynamics in d > 2. In the present work we prove 
rigorously that Equation (0) is a necessary and sufficient condition for the existence of 
product measure steady states for ZRPs, and as a consequence that the system described 
by (jSJ) has no product measure steady states. This condition in turn is a special case 
of a condition on a class of systems which we call Generalized Zero Range Processes 
(GZRP), in which we also allow transitions in which more than one particle moves 
at a time, although we will assume that the number of particles movimg in a single 
transition is bounded. The rate now depends on the number of particles v which move 
in the transition as well as the number of particles n at the site before the transition, 
and so the rates are given by a function a q (u, n) with some v max such that a q (n, v) = 
whenever v > v max . The classical ZRPs discussed above are a special case with v max = 1 
and a q (l,n) = a g (n). We prove that a necessary and sufficient condition for the GZRP 
to have product measure steady states is 



for some non-negative g q (y) and f(n) with f(n) < oo. This has a clear similarity to 
(J7J), and when v max = 1 © reduces to ((TJ). 

We will prove Equation (JSJ) in the course of finding a weaker result for the 
continuous-time Mass Transport Process generalized to dimension d > 1 and to 




(7) 



a±i(n) = a[l - 6 nfi ] 




(8) 




(9) 
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transitions in all directions. We show that these systems have product measure steady 
states when a q (fi,m) = g q (fi)p(m — n)/p{m). This condition is also necessary under 
certain conditions (generalizing (|ZJ) to higher dimension) and we conjecture that this is 
so in all cases. 



2. Factorizability in the Mass Transport Process 

Let Pa(zzM) be the time-dependent probability density of finding the system in a 
particular configuration m with mass m; at site i G A, m; G (0, oo). The master 
equation describing the evolution of Pa(zzM) is 

dP A (rn,t) 
dt 



^ ( ~y] I dfjLa q (ji,mi)P(rn,t) (10) 



dfxa-qifx, m l+g + jj)P(m i,? '", t) 



where 



m- q ^ = { nij — fi, j = i (11) 
rrij + //, i = i + q 

A stationary state of the system is a distribution Pa(zh) such that dP\(m, t)/dt = 
whenever Piv{Ek,t) = P\ (m) , or equivalently 

rrrii _ nmi 

^2^2 di2a q (/j,,mi)P A (m) = ^^2 d/2a^ q (iJ, } m i+q + /j,)P A (m t ' q ^) (12) 

ieA q ^° ieA q ^° 

We wish to find conditions under which there is a Pa which is factorizable, that is 
which takes the form of 

PA(ni) = l[p(m i ) (13) 

ieA 

Assuming that a q (n,m) > for all m > and < /x < m for at least one q of 
each pair of opposite directions, the system in a finite torus A will have a unique steady 
state corresponding to each value of the total mass M = m i- Any linear combination 
of such states is also a solution of (fT2"|) . Given a factorizable steady state, states of 
definite total mass can be obtained by projecting P onto the set of configurations with 
a particular value of M in analogy with Equation (J2J) . 

Assume that there is a factorizable steady state as in (|13j) . Let p(s) be the Laplace 
transform of p(m), and let 

/■oo 

(f> q (li, s) = [l/p(s)] / dme~ srn a q (fi,m + fi)p(m + ft) (14) 
Jo 

Note that, since a q (fj,, m) = for m < /i, 

/■oo 

/ dme~ sm a q (fi,m)p(m) 



Product Measure Steady States of Generalized Zero Range Processes 5 

poo 

I dme~ sm a q (p } m + p)p(m + p) = e~ s ^4> q (p } s)p(s) (15) 
Jo 



We also have 

) pm 

dme~ sm / dpa q (p,m)p(m) = I dp I dme~ srn a q (p,m)p(m) 



POO POO 

= / dpe~ Sfl / dme~ sm a q (p,m + p)p(m + p) 
Jo Jo 

POO 

= / dpe- s ^ q (p,s)p(s) (16) 
Jo 

Multiplying both sides of (|T2*j) by Yli e~ Simi and integrating over all mj, we obtain 

(\ poo pim 
Y[l'( s j) / dnH dpa q (p,mi)p(mi)e- s ' m ' 

t\ poo poo poo 

II / dmi / rfm *+« / d A* 

x a- q (ji, m l+q + /i)p(mi - p)p{m i+q + ^ e -^- Sl+q rn l+q ^ 

Rewriting fTTj) with the aid of (|T5|) and (fT6j) and cancelling common factors, we 
obtain 

POO POO 

dfi<f> q (ji, s^e-"^ = / dp<j)_ q (p, s i+q )e- s ^ (18) 

Equation ()18j) will be satisfied if (though not only if) 

<f>q(P, a) = g q (p) (19) 

In this case Equation (JT3jl gives 



oo 



dte sm a q (p,m + p)p(m + p) = g q (p) / dme sm p{m) (20) 

Jo 

which by uniqueness of the Laplace transform gives 

a q {p, m) = g q (p) P m ^ (21) 

Equation (|2ip is a generalization of the comparable formula for the unidirectional 
case jZj. In this case and in all other cases where, for each q £ E, either a q = or 
a_ 9 = and hence either <p q = or 0_ g = 0, there is in Equation (|18|) only one term 
which depends on each pair m», m i+q , and in order for the equation to be satisfied it 
must depend on only one of them. This happens only if (|19j) holds for that q, so in these 
cases Equation ()21|) gives the only possible rates for which there is an invariant product 
measure. 

Although in general Equation (|T9*|) is not the only way of satisfying Equation ([18)1 . 
solutions of this equation only correspond to realizable dynamics when p and a q are 
non-negative and normalizable; the resulting restrictions on (p q from Equation ()14j) are 
such that it seems unlikely that there are reasonable (indeed any) rates, other than 
those in (|21|) . which satisfy all of these conditions. 
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Dynamics for which the system has a factorizable steady state can be found 
by beginning with some suitable (positive and normalizable) p(m) and then defining 
a q (p,m) via (|21j) . For example let 

p c (m) = ce- cm 9{m) (22) 

where 6 is the Heaviside step function. The possible transition rates corresponding to 
P(m) = YlPc( m ) are 

a q (m, fx) = g q (fi)e Cf *9(m - //) = g q {p)6(m - //) (23) 

where g q are arbitrary non-negative integrable functions, i.e. the rates a q (p,m) are 
independent of m as long as p < m. 

3. Reverse processes 

The existence of product measure states of the Mass Transport Process is related to the 
nature of the reverse process of the Markov process defined by ([lip. The reverse of a 
Markov process defined by transition rates K(m — *> rnf) with stationary state P(m) is 
defined by rates K* (mf — > m) given by 

, , Kim — > m')P(m) , , 

fT(m' -> m) = v — . — ; y —> (24) 

P(m') 

Examining ()24|) . it is clear that the reverse process of the MTP involves only 
transitions in which a mass fi moves from a site i to a site i + q, q £ E. Denoting the 
rates of these transitions by a* q (rn, pi), and assuming that P(m) is a product measure 
with single-site weights p(m), (|2~3j) becomes 

a g (/i, + n)p{mi + n)p(m i+q - /i) = a* +g _ g (/i, m)p(m)p(m i+q ) (25) 

from which we can see that a* depends only on q, mi, and m i+q . We further note that 
a* is independent of m i+g , and so defines an MTP, if and only if a q satisfies Equation 
(|21|). We therefore find that the statement that Equation (|2~T|) is a necessary condition 
for an MTP to have product measure steady states is implied by the statement that the 
reverse process of any MTP is also an MTP. 

4. Factorizability in Generalized Zero Range Processes 

With mass at each site restricted to an integer particle number n^, we can reproduce 
the analysis in the previous section up to Equation (fTHjl . Denoting the vector of 
occupation numbers by n, and the transition rates by ct q [y, n), the stationarity condition 
corresponding to Equation (|T2*|) is 

rti 

^)^a(^) + + ^a(z^)) = (26) 

ieA q<=E u=l 
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Suppose P is factorizable, 

h(n) = Y[p(rn) (27) 

where p(n) is the probability of having n particles at a given site. Then define the 
generating function (discrete Laplace transform) 

oo 

p(*) = 5Vp(n) (28) 

n=0 

and let 

M», z) = ^o^ q ^n + u) P (n + ,) (2g) 

p{z) 

Note that (j) g (v,z) > for all is, z > 0. The counterpart of Equation (fTS|) is then 

oo oo 

E ^2 z iM i/ > z i)= E E ^ -^+9) ( 3 °) 

We now exploit the assumption that transitions occur only for v < v max . Then 
choosing some j G A and q G E and taking the fcth derivative of the above expression 
with respect to Zi and Zi +q gives 

E T^ZTV } & z ^ + E TTTZTv %£*T *) = ( 31 ) 

For = v m ax, we have 

^^{Vmax , 2j) + 0^ (*W , Zi+q) = (32) 

For (J32)) to hold for all z^ and Zj+q, both terms on the left-hand-side must be constant, 
and thus the functions 4 > ± q (is max , •) are polynomials of degree v max \ being non-negative 
they must have non-negative leading terms. Equation (J32J) states that pairs of leading 
terms of these polynomials must add up to zero and so each must be zero, and therefore 
the functions (f) q {y maxi ■) are polynomials of degree at most v max — 1. 

Now setting k = v max — \ we find by the same reasoning that the functions 4> q {y max , •) 
and 4> q {i>maxi ■) are polynomials of degree at most v max — 2. Proceeding in this manner 
we find 

<t>q(v,z)= 9q(y) ( 33 ) 
as a necessary as well as a sufficient condition. Referring to the definition of 0, this 
implies that 

a q {u, n) = g q {v) -— (34) 

p[n) 

is a necessary and sufficient condition for the existence of a product measure. 
In the case where v max = 1 and a q (l,n) = a q (n), this condition becomes 

a q (n) = c q -— - (35) 

p[n) 

This is what we referred to above as the classical ZRP, with the well-known stationary 
measure [21 E] discussed in the introduction. 
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5. GZRPs on infinite lattices 



It remains to be established that there are not GZRPs on infinite lattices which have 
product measure stationary states which do not correspond to stationary states of the 
corresponding GZRP on a finite periodic lattice. Let P(n) be a product measure with 
single-site distribution p(n) which is stationary for rate functions a on Z d , and let A be 
a finite box in Z d which is at least two sites thick in all dimensions. Denote by n A the 
configuration of the system inside of A, and let P(n A ) be the marginal distribution of 
this configuration. Then we have 

d ni Ui+q 

p(n A ) = - Y Y Y n i) p ^h) + Y Y Y a i( u > n * + v ) p h± 

ieA q€E u=l ieA qeE u=l 



dt 



oo oo 



~Y Y YY a o( u > n ") p (-^ p ( n ") 

i£dA q&E:i+q£A v=X n=v 
oo n 

+ E Y YY a ^ p< ^ 9 np(n) 

i£dAq£E:i+q£An=l u=l 

= (36) 
where OA = {i G Z d \ A\(3q G E)(i + q G A)} and 

{n k , k {i,i + q}HA 

n k + v, k = ieA (37) 
rife — v, k = i + q G A 

By repeating the procedure used above with Equation (J3TI|) in place of Equation (ffilj) . it 
can easily be seen that a and p must satisfy Equation (J3T|) and so p gives a stationary 
product measure for the rates a on any finite lattice. 



6. Conclusion 



We have shown that there is a straightforward necessary and sufficient condition, 
Equation Q, for a generalized Zero Range Process to have a product measure steady 
state. For Mass Transport Processes, we have found a condition, Equation (fTSj). for the 
existence of a product measure steady state, which is considerably more opaque than in 
the GZRP; it is not clear that this is equivalent to the sufficient condition expressed in 
Equation (J2T]) . the counterpart of the condition we have obtained for GZRPs. We have, 
however, presented some reasons to believe that it is. 
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